Seminal magnetic fields from Inflato-electromagnetic Inflation 
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We extend some previous attempts to explain the origin and evolution of primordial magnetic 
fields during inflation induced from a 5D vacuum. We show that the usual quantum fluctuations of 
a generalized 5D electromagnetic field cannot provide us with the desired magnetic seeds. We show 
that special fields without propagation on the extra non-compact dimension are needed to arrive 
to appreciable magnetic strengths. We also identify a new magnetic tensor field Bij in this kind of 
extra dimensional theories. Our results are in very good agreement with observational requirements, 
| in particular from TeV Blazars and CMB radiation limits we obtain that primordial cosmological 

magnetic fields should be close to scale invariance. 
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& . I. INTRODUCTION 

on 

Magnetic fields seem to be ubiquitous in the universe. Observations have well established the widespread presence 
of magnetic fields in the universe (Tj-Q- Cosmic magnetism has been verified of /zG-order strength in galaxy clusters 
j— |. and also in high redshift protogalactic structures. Recently, in particular, Kronberg et al and Bernet et al reported 
organized, strong B-fields in galaxies with redshifts close to 1.3 [H All these seem to suggest that magnetic fields 
similar to that of the Milky Way are common in remote, high-redshift galaxies. This could imply that the time needed 
by the galactic dynamo to build up a coherent _B-field is considerably less than what is usually anticipated. On the 
other hand, the widespread presence of magnetic fields at high redshifts may simply mean that they are cosmological 
(pre-recombination) in origin. Although it is still too early to reach a conclusion, the idea of primordial magnetism 
t^J" \ gains ground, as more fields of micro-Gauss strength are detected in remote proto-galaxies. Recent reviews about 
J> . primordial magnetogenesis can be found at [Toj . 

Our work focuses in studying the production of primordial magnetic seeds during an inflationary period, which in 
turn is induced by the immersion of a particular 4D hypersurface, given by a de Sitter spacetime, in a 5D vacuum 
space defined on a 5D Riemann-flat extended de Sitter metric. The theory that motivates such scenario is the Induced 
Matter Theory (IMT) [ill [l2j . In this theory inflation can be recovered from the Campbell-Magaard theorem[IM3, 
which serves as a ladder to go between manifolds whose dimensionality differs by one. This theorem, which is valid in 
any number of dimensions, implies that every solution of the 4D Einstein equations with arbitrary energy momentum 
tensor can be embedded, at least locally, in a solution of the 5D Einstein field equations in vacuum. Because of 
this, the stress-energy may be a 4D manifestation of the embedding geometry. Physically, the background metric 
there employed describes a 5D extension of an usual de Sitter spacetime, which is the 4D spacetime that describes an 
inflationary expansion. 

We will study the simplest inflationary model, that gives us a de Sitter epoch. Electromagnetic studies in this 
context have been developed previously by [18-20]. To introduce the electromagnetic effects we consider a massless 
vector field A a of five components. The component normal to the hypersurface has properties similar to a scalar field, 
so its spectrum, with the Coulomb Gauge, we will see it can only be scale invariant. 

In the present paper we consider two types of fields. The first with propagation outside the hypersurface, and the 
second with the propagation confined to the hypersurface. This last fields need to be generated by some brane source 
since they produce a discontinuity in the stress tensor. Such distinction is very important since the fields that can 
propagate outside the hypersurface are strongly suppressed than the usual 4D photon, giving rise to a bluer spectrum. 
The other fields, in contrast, can give invariant spectrum for magnetic fields. Our results in this way are very similar 
to 0. 

Additionally, we analyze the electromagnetic fields measured by 4D observers in the effective 4D hypersurface. 
When we extend our definitions of the relativistic electric and magnetic fields to more dimensions, we obtain new 
gravitoelectromagnetic fields of scalar and tensorial nature. In our knowledge this is the first time such components 
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are identified in an electromagnetic theory of Brane Worlds or Induced Matter Theory. Finally, it is well known that 
the spatially flat FRW universe is conformal flat and the Maxwell theory is conformal invariant, so that magnetic 
fields generated during inflation would come vanishingly small at the end of the inflationary epoch. However, in 
spatially open FRW models magnetic fields are sufficiently strong to seed the galactic dynamos. The reason is that 
conformal flatness of the 4D background metric is a global property in spatially flat FRW spacetimes, but a local one 
in spatially open or closed geometries {22H24I I . 

This is a very important problem of inflation to explain seminal magnetic fields. The possibility to solve this 
problem relies in produce non-trivial magnetic fields in which conformal invariance to be broken. The conformal 
invariance of Maxwells equations in four dimensions can also be broken if an embedding into a higher dimensional 
space-time with time- varying extra spatial dimensions is considered. In particular, in Brane Worlds or Induced Matter 
theory conformal invariance is naturally broken, which make possible the super adiabatic amplification of electric and 
magnetic field modes during the early inflationary epoch of the universe on cosmological scales. 



II. ELECTRODYNAMICS A RIEMANN-FLAT EXTENDED DE SITTER SPACETIME 



We begin considering the action of a free abelian gauge vector field in five dimensions 



S = -±Jd 5 xV^~gF CD F CD ; (1) 

where Fqd = VcAd — V^Ac = 9cAd — doAc is the antisymmetric field tensor. The Maxwell equations of motion 
in five dimensions are 

V C F CD = 0. (2) 

We introduce these fields in a background 5D spacetime in vacuum state. Using some ideas of the Induced Matter 
Theories we want to obtain a hypersurface that undergoes an effective inflationary period. In this sense we do a the 
semiclassical expansion of the vector fields 

A C = A C + 5A C , (3) 

where the overbar symbolizes the 3D spatially homogeneous background field consistent with the fixed background 
homogeneous metric and SA C describes the fluctuations with respect to A c . A more detailed explanation can be 
found in [28[. This space will be a 5D extension of a 4D de Sitter spacetime, which is very relevant to inflationary 
cosmology, and is given by a 5D (canonical Riemann-flat) metric with a line element [2^ 

ds 2 = 4' 2 dN 2 - ^ 2 e 2N dr 2 - d^ 2 , (4) 

where N is a dimensionless time like coordinate related with the number of e-folds during a de Sitter inflationary 
expansion. The 3D space like cartesian dimensionless coordinates dr 2 = dx\ + dx\ + dx^ when taking a constant 
foliation coincide with the 3D space coordinates of the effective 4D hypersurface. The extra (non compact) space-like 
coordinate ip has spatial dimension. Instead of these coordinates, we shall use new conformal ones [see appendix l[A"|)]. 
such that 

ds 2 = b(w) 2 [a{r/) 2 (dr] 2 - dR 2 ) - dw 2 ] , (5) 

where the scale factors b(w) — e H ° w y a(rj) = — -gi— are dimensionless. The components of the coordinates (77, R, w) 
have spatial dimensions. The foliation ip = tpg corresponds with w — Wq = in these coordinates, such that b(w a ) = 1. 
The volume factor of the manifold is \A 5 )o = b 5 a 4 = e r . 4 4 , and the non zero connections are 

r^-rr 1 ^, rl B = e- 2H ° w H g AB , r^ = H 6^. (6) 



1 In our conventions indices "A,B,C,..,H" run from to 4, Greek indices run from to 3 and latin indices "i,j,k,..." run from 1 to 3. 
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The equations ([2]) in these coordinates are 
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Hq drj \dw 
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(7) 
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Hq \dw 2 dw 



A, + rfd&Ai + -^dj 



d_ 

dw 



H )A 4 = 0, 



(8) 



07] z rj or] J dw dw \r) or) t 



(9) 



where d 2 is the ordinary Laplacian operator. 

We adopt the usual Coulomb Gauge for the study of the electromagnetic effects on the hypersurface: A$ — y 
diA 1 = 0. The equation for the temporal component Aq remains as a constraint equation for A 4l 



d_ ( d_ 

dr] \dw 



H )A 4 = 0. 



(10) 



Notice that in a 5D Minkowskian metric: dS 2 — r\ a \ > dx a dx b , the connections vanish and the field equations remain 
decoupled after the gauge choice. In this case the universe does not expands so that the inflaton field A 4 = 0. 
To continue we propose the ansatz 



A^{r 1 ,K.,w)=e mH ° w A m {r,,K), 



(11) 



where to can be any complex number. In particular, to make a Fourier expansion base functions, we shall consider 
the allowed values to = — 1/2 + ip, where p = P/Hq are the wave numbers of the field in the w direction per 
unity of Hq. The reason to make it relies in that y/g = e 5 #o«>^-4 an( j FabF ab oc e~ 4HoW , so that we obtain 
J dr] J d 3 x dw e H ° w (- ■ ■ ) in the action. This exponent should be suppressed in order to address a correct Fourier 
expansion in terms of e lP '"+ lk - R . 

It is also important to consider particular solutions of perturbative order, without the oscillatory regime with respect 
to the extra coordinate. These solutions are described by the real values of m. In any case the equations are 



(1 + m)—A 4m 



0, 



__^]-m(m + l) 



Ah 



3 m + (m+l)d 3 A 4m = 0, 

d 2 2d 
drj 2 rj dr] 



- d 2 A im = 0. 



(12) 
(13) 
(14) 



We notice from (|T2"j) that the only possible mode for A 4m that can exist is to = —1. The others yield d 2 A 4 = 0, such 
when we suppose the absence of non trivial 3D surface conditions. For m = — 1 the equations of motion 



that Ai m 
are 



drj 2 



d 2 \A 



3 m— — \ 



= 0, 



d 2 2 d „,\ „ 
— —-d 2 )A 4m= _ 1 = 0. 
dr] rj dr] J 

Furthermore, the solutions with to 7^ — 1 remain decoupled since Aj m j_i 

d 2 

— r - d 2 ) - m(m + 1) Ai 



0, 
0. 



(15) 
(16) 

(17) 



2 The unique nonzero background field is the inflaton which drives inflation and has a constant expectation value on the effective 4D de 
Sitter hypersurface [2^|. Hence, in the following we shall denote the fluctuations A c = 5A C 
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III. ELECTRIC AND MAGNETIC FIELDS IN 5D 

We need to define the electric and magnetic fields for the observers that belong to the hypersurfacc. 

When we extend the electromagnetic theory to five dimensions the number of independent components of the 
antisymmetric electromagnetic tensor Fqd is 5(5 — l)/2 = 10. But the interesting fact is that the dual tensor is of 
third rank. 

J~cde — -^]cdeabF ab ', (18) 

where the space volume tensor tjcdeab — \fg ccdeab is of 5th rank. If we want to define the electric and magnetic 
components seen by observers, we can apply the same definitions as in 4D. In this way we obtain a 5D vector electric 
field 

E c = u D F CD , (19) 

and a 2-tensor antisymmetric magnetic field 

B C D = \riCDEABU B F EA . (20) 

In particular, for a 5D comoving observer defined by u A = a _1 6 _1 (l, 0, 0, 0, 0) [in conformal coordinates ((SJ)], we 
obtain the following expression for the components of this antisymmetric magnetic tensor 

B( DD ) = 0, (21) 
Bod = 0, (22) 

Bu = ^y/ge 4ijko u°F^, (23) 

Bij = Vgeijk40U°F kA , (24) 

where we have distinguish the temporal, the 3-spatial and the extra coordinate of the hypersurface. The number of 
independent components is six: 3 belong to B4& playing the role of the usual 3-magnetic vector field; while the other 
3 components belong to a new 3D antisymmetric magnetic tensor. 
The electric field is also decomposed for this observer as 

£o = 0, Ei = \d Ai, E i = \d A i , (25) 

a o ab 

there are 3 components Ei with the same properties as the usual electric field, and an extra component E± that we 
call the scalar electric field. Then the 10 components of the tensor F a b decompose into: 3 of the usual electric vector 
field, 3 of the usual magnetic vector field, 3 of a new 3D magnetic tensor field By and 1 of a new electric scalar field 
Ei. 

A. Quantization in the Coulomb Gauge 

Once we apply the Coulomb Gauge the dynamical fields are Aj y Aj. To perform a canonical quantization we 
should impose commutation relations between the field and its conjugate canonical momenta. From the 5D action 
((TJ) we define the conjugate momenta n B = = \fgg 0Q g BC Ac In the coordinates (77, R, to), they take the form 

= &H 3 a(r/)V^,o, ^ = b(w)a( V ) 2 A 4fi . (26) 
Therefore, the commuting rules at equal times are 

[A 4 (r/,R,u;),7r4(?7,R',u/)] = b 3 a 2 

dA n 



A 4 (77,R,w),^(r/,R,u/ 



A < fa,R,u;),^(» J ,R> / ) 



= id^{R-R')S(w-w'), (27) 
= iS(w-w r )6 { l r AR-Il'), (28) 



[A j (rj, R, w) , TTi (77, R', w')] = b 3 a 2 

where 

is the 3D transversal Dirac Delta distribution, which takes into account only wavenumbers which are perpendicular 
to the direction of propagation in order to obtain a correct quantization in the Coulomb Gauge. 
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B. 5D Fourier decomposition 

The 5-vector A c has five polarization states. Considering solutions from ([2]) that propagate in all directions 

Ac(s,k,p\r,,R,w) = e Cs (k,p)A S}k , p (r))e ik -' R - + ^ HoW , (29) 

where s labels the polarization states, k is the comoving wavenumber in the 3D-spatial space of the hypersurface and 
pHo is the wavenumber in the extra direction w. To accomplish with the Coulomb Gauge requirement we introduce 
the next orthonormal base of polarization vectors: 



b^a' 1 (1,0,0), 



0,-^,0), < 



& _1 (0,0,1) 



(30) 



Here, A = 1, 2 represents the transversal polarizations to the 3D-spatial propagation of the wave, they yield e\ ■ e\ = 1 
(no summing over A). With this choice the transversal vectors yield = £*a • k = 0, that it is the Coulomb Gauge 

in momentum k space . The factor b~ l a~ l is important to define a vectorial base. The completeness relation derived 



5>-\ 



A B 



AB 
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(31) 



s=0 



where the factor h ss is matrix with the components of the 5D Minkowski tensor metric. The previous expression 
projected to the 3D space is 



k l k l 



£e A (k)e jA (k)+<^ 



A=l 



(32) 



This expression is very important to build the commutation rules. Finally, the projection in the extra coordinate 
simply yields 

h 44 4 =4 e* =4 = g 44 = -e- 2HoW . (33) 



Using the last considerations, we expand the field Ai(rj, R, w) in the solutions (129 

//.. :<■ I ti' I: tin 

Aj(rj, R, w) = e 



dp 



(27r)3/2 (27r) i/ 2 E^^f) [bx(K P )e^^A k , p (v) 
+bl(K P )e^ H ° w -*- R Al p (i 1 ) . 



(34) 



The additional factor e - H o w / 2 comes from the expansion in functions e <p - ff °' u ' +tk ' R in such a way that the integrand 
remains normalized. 

Because of the factor b~ 1 a~ 1 in the polarizations, the function that yields the equation of motion (|17[) is: Ak, P {ii) = 
b(w)a(r/)Akp(r]). Then the mode equation is 



d Akp 
drf- 



„2l 



A kp = 0. 



(35) 



Notice that the value m = — 1/2 + ip yields — m(m + 1) = |m| 2 = 1/4 + p 2 for the solutions that propagate on the 
extra coordinate. The creation and annihilation operators satisfy the usual relations 



b x (k,p),b[,(k',p')] = 6W(k~k')6(pHo-p'Ho)5 X y, 



[6 A (k,p),6 v (kV)] = k(k, P ),&t,(k>') 



0. 



(36) 
(37) 



Finally, the normalization condition for the temporal modes Ak P (r]) that comes from the commutation relations is 



a 2 b 2 W[A kp {r])} = i, 



(38) 



where W[f] — //'* — /*/' is the Wronskian and (') denotes the partial derivative with respect to the conformal time 

TJ. 

We have seen that the only scalar mode that propagates is Ai m= _i, this means that an analog quantization in 
functions e lpff °"' +lk ' R won't work in this gauge for A4. In its place we suppose small perturbative fields that only 
propagate in the 4D hypersurface, decaying exponentially outside. The origin of such fields is not addressed in the 
present paper. 
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C. Vector modes 



In the following we will distinguish between modes that propagate in all space directions, including the normal to 
the hypersurface w, where m = — 1/2 + ip, and modes that only behave like plane waves in the hypersurface, with 
met 



1. Extra dimensional propagation of the modes 
The solution for ([55)1 of the modes that propagate in all directions is 

A kp (v) = d HkiP) (kT,)^ 2 n^(kT,) + d2 {Kp) {kT ] ) 1 ' 2 nf\kr 1 ), 



(39) 



where i-L p 1 \krf) = (^J p + iY p ^j {krj) and H p 2 \kr]) = (^j p — iY p ^j (krj) are the Hankel functions with imaginary order 

v = ip, J p (x) and Y p {x) are the first and second kind Bessel functions with imaginary order defined real in the 
following way [30( [the reader can see some properties of the Bessel functions with imaginary order in the appendix 



J p (x) 



cosh 



r^Re [J ip (x)} , 



Y p (x) = 



cosh 



-^-Re [ljp(x)] 



(40) 
(41) 



Here, we take into account the real part of the usual decomposition of the Bessel functions. The quadratic amplitude 
of these modes in the large scales is 



Ak P A* kp \iR ~ - 
p 



tanh 



(?1 
\ 2 



COS" 



pin I ^ 



7 P 



pin 



7p 



tanh(^) 



(42) 



In the figure ((5} we show the amplitude of the modes in the effective 4D hypersurface versus p for two different 
wavenumbers k, during the inflationary epoch. The amplitude is not divergent for p — ¥ 0, since in this limit case we 
have (tanhp) p _i.o — » p, (coshp) p _i.o —> 1 and (sinhp) p ^ — > Pi so that 



hi) 



In 



(43) 



where 7 ~ 0.5772.. is the Euler constant. On the other hand, when p ^> 2/n one obtains (tannp^^/^ — ^ 1> while 
the quadratic sines and cosines have the same argument reducing to unity. Finally, we obtain 



• / 4fcp>2/7r-^fep>2/ 7 r \IR 



rj a 



P pHo 



(44) 



Other important quantity is the amplitude of temporal derivatives of the modes. In the long wavelength limit it takes 
the form 



•A'kp<Akp\iR 



ip 
V 



tanh 



2 J 



sin 



pin 



cos 



pln(M) 



7 P 



tanh (SE) 



(45) 



In the figure (J6j> we show the amplitude of the derivatives of the modes for two different wavenumbers k. When p = 0, 
we obtain 



•Afcn-/4.fr! 



8H Q 



while for p ^> 2/tt, we get 



feO-^feOU-R 



•A'k P <Ak P \iR 



TTTj 



4p 



ApH a{rj), 



(46) 



(47) 



so that in both limits grows as Hq a(rj). 
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IV. ELECTRODYNAMICS ON A 4D HYPERSURFACE 



With the aim to study the effective 4D dynamics of the fields without propagation in a de Sitter spacetime, we shall 
consider a static foliation w = Wq on the metric ([S]). This hypersurface is relevant to describe a de Sitter expansion 
of the universe in the early inflationary epoch. 



A. Effective 4D quantization 



Once we make the foliation on the fifth coordinate w — wq and the modes are quantized [see appendix (|B])]. we 
obtain the modes Ak m (j]) = b(w) a(rj)Ak m (j]) that yield the equations of motion (|17|) 



d 2 A k 



dr\ 2 



, 2 to(to + 1) 



A 



o. 



(48) 



which comply with the normalization condition 



b 2 a 2 W[A km {ri)] = i, 



(49) 



that is compatible with the commutation relationships (|B6 



B. Infiaton modes on 4D 



Once fixed the Coulomb Gauge, and with our ansatz there is only one possible solution for the extra component 
A4: the one with the value m = — 1. We call this solution as <j>k, that behaves as a massless scalar field in 4D. Using 
the equation (|4"g]) we obtain the general solution 



(50) 



From the condition (j49|) we can set one constant to zero C2 = 0, so that ci(fc) 



In the infrared limit we obtain 



°k\ IR 



Ho 
Tib 



4 ,a 3 .3,_ a 

-k*rf + t-k 2 



(51) 



This spectrum is exactly scale invariant, something we would expect because we are dealing with 4D hypersurfaces 
w = that suffer a de Sitter expansion. In this case the amplitude of the modes, \4>k\ = yj<j>k<j>k, is constant 



3ff 
2nb' 



(52) 



in agreement with what one expects in a de Sitter expansion for a massless scalar field. In this sense this special 
solution has the appropriate spectrum of the infiaton fluctuations. 



C. Vector modes confined to the 4D hypersurface 



Lets turn to study fields where their modes are confined to propagate in the effective 4D hypersurface, which de- 
scribes a de Sitter expansion for the universe. We proposed an ansatz for such fields as A c (i], R, w) — e mHow A^ l (rj, R), 
depending exponentially in a parameter m6l. This choice, however, can be problematic when the parameter m < 0. 
In this case the fluctuations blow when w — > -co. In the same way, if m > the fluctuations blow as w — > 00. To 
avoid this problem we can suppose that these fluctuations reach a maximum in the 4D hypersurface, and therefore 
decay exponentially at both sides. A possible symmetry to choose is Z2, common in brane theories. In this sense, for 
example, if m < 0, we have the solution e ~ mH o w f or w < 0. However, the solution will be e mH ° w for w > 0. In this 
way the function is continuous in w = 0, but not its derivative that suffers a discontinuity 2mHo- Apart from this, 
it appears a complication since the equations of motion (1481) change when we go from one side to the other of the 
hypersurface. This is due to the fact they are not invariant under the transformation m O —m. A closer inspection 
shows us that they are invariant under the change m ■<-> —(to +1). This means that the m-mode is equivalent to 
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FIG. 1: Dependence of the modes outside the effective 4D hypersurface using the shift symmetry m o — (m + 1) 



the — (to + l)-mode. Then, we use the solution e -(" l + 1 )#o"' wnen m < _i f or the interval w < and e mH ° w for the 
interval w > 0. On the other hand, when m > the solutions should be e mHow for io < and e -( m + l ) H ow f or w > 0. 
This idea is sketched in the figure (JXJ) . In the interval —1 < to < we don't have this symmetry and the solutions 
would irremediable blow. For this reason they remain excluded. 

The solution written in terms of Bessel functions, reads 

J-i- m (kv)\- (53) 



The quantum limit of this modes is 



A km \uv — > -F=e lkr >, krj — >-oo. (54) 

V 2k 



Using the normalization condition and the asymptotic form of the Bessel functions we find 

ci(fc, m) - -c 2 (fc, m) e- l7r (^ +m ), (55) 

, - (1— m) 

C2(fe,m) = W77 ; v-. (56) 

y ' ' V 4fc cos(to7t) 
Once defined the modes for all scales we can find their large scale limit for (krj — > 0) 

AmWUil - Y W1 ' -— k-*(kr,)- m + (57) 

22-™ r (J — m) cos(m7r) 



/7T e l7r 2 



2 |+m p (I + m ) cos((m + 1)^ 



k-?{krj) m+ \ (58) 



This is the long wavelength solution for the modes to < — 1 when u> > 0. It is also the solution for to > when 
it) < 0. It is interesting that the solutions do not change when we transform — m (m + 1), in this case the first term 
shifts to the second and the second to the first. It is useful to introduce the following dimensionless function [2l[ to 
calculate the amplitude of the modes. 

^ = 2 2M+l r 2(l +A1 )cOS 2 (^)' (M) 



where 



-to, to > — — /i = to + 1, m < — — . (60) 
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F(-3) = 4.5 



-4 -3-2-1—0 1 2 3 

2 

FIG. 2: The amplitude J-(m) for modes without propagation outside the hypersurface. The invariant spectrum for magnetic 
fields will be achieved for m = —3. The dashed sector is forbidden so the fields are fluctuations that do not diverge in the extra 
coordinates infinities. 



The function only diverges in fi = —1/2, dividing the function in two well behaved branches. Previously we noted 
that only the modes with m < — 1 and m > have the symmetry to avoid divergences and infinities. This means that 
for this modes the function is well behaved. 

The amplitude of the modes in the hypersurface w = is 



Another useful quantity is the temporal derivative of the modes Ak m - We find from (|57[) 



A' km (v) - Vnk* 



23-T (i _ TO ) cos (7rm) 



(krj) 



+ 



+ to) cos(7t(to + 1)) 



(k V y 



Again we define a dimensionless function to describe the amplitude of the modes 

7T(0 + 1) 2 



where 9 depends of m in the following way 



2 3 + 2e r 2 (f + 6) cos 2 (t:(9 + 1))' 



-(m + 1), m > 



= m, m < — - 



(61) 



(62) 



(63) 



(64) 



The amplitude reads 



A' ke A' k * e ~ kG{6)(kri) 



21) 



(65) 



V. ELECTROMAGNETIC ENERGY DENSITIES 



The stress tensor we find from the action is 

Tab = -9 CD F AC F DB - \g A B9 CD g EF F DF F CE . (66) 

In the same way as in the usual 4D Maxwell theory, where Too = —(E 2 + B 2 ), we can introduce an expression for the 
temporal-temporal component of the stress tensor as the sum of the individual energies of each physical field 

Too = T E j + T B j + + . (67) 
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Lets write explicitly each of the previous contributions to the energy density as seen by a 4D comoving observer. The 
vector magnetic part of it, is 

T^>oo = -fVV^. (68) 
Since BnB Al = g iA BiB l and Bi = eijka~ 1 djAk, we get 

B it B 41 = b- 4 a- 4 [djAkdjAk - djAkdkAj] , (69) 

so that 

T^oo = -^B,^ = ia 2 B^. (70) 

This is the same result that obtained with the usual 4D electrodynamics. The contribution due to the electric field 
(vector) part is 



and the electric (scalar) part 



T^ 00 = ^-E 4 E\ (72) 



Finally, the contribution due to the tensor magnetic part is 



T^oo = -^-B^K (73) 



The energy density is defined as p = — (T§) . In our coordinates it takes the form 

P=\ - E 4 E 4 + B it B il + Bu&t) . (74) 

All the terms are positive, EiE % = —b~ 2 a~ 2 Ef and E4E 4 = —b~ 2 E\ [in comoving coordinates {5}]. The previous 
result in any coordinate system reads 

9 = g (-EaE a + B CD B CD ) . (75) 

A. Energy density due to Magnetic vector Bi 

We can calculate the contribution of the magnetic vector to energy density defined as Piba = — {T^ Bi \)' 

(T^° ) = -l(B ii B ii ). (76) 
Using the expansion in the modes A^ p , we obtain its derivative in the direction j\ 

d n A n = f ,f* * fej^ R +«A fep + HkjJbxe -W^o W ) Al \ . (77) 

J (2tt)2 (2tt)2 ^ V / 

Its vacuum expectation value is given by 

/ f ]3 l. fj n J 2 
—r \A kp \ 2 Y, k n *fc ex h e Xh . (78) 
(27r)2 (27r)2 r-~ 



A=l 
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So that (T^^q) can be written as 



(T 



1 



d 3 k dp 
2& 5 a 4 J (2tt) 3 2^ 

2 

x S ' (kjkjei\£i\ ki£i\kj€j\) , 

A=l 



(79) 



Using the completeness relation projected to the 3D space and the transversal condition k^x — 0, the second term 
in the previous expression vanishes, while the first one holds 



2 

E 

A=l 



(kjkjCixeix - kieixkjCjx) = 2k 2 b 2 a 2 



(80) 



Since the solution depends on the absolute value of p, we can write the integral as dp = 2 dp, so that the 
magnetic contribution to the energy density is 



P(Bi 



1 



2tt 2 6 3 
1 

2ir 2 b 5 



A 



kp 



dk dp , - 
T — k 

k 7T 



dk dp 5 \Akpl" 2 
k 7r a 4 



The energy density stored at a certain scale fc for unity of p is 



dkdp 



P{B t ){v^.p) 



1 



2tt 3 6 5 a 4 



\A kp \ 



Recalling ((32) we arrive to 
d 



dkdp 



P(Bi){v,k,p) 



1 



k 4 



2tt 3 pH b 5 a 5 



tanh 



(f) 



cos 



pin ( — 
1 2 



1p 



phi 



1p 



tanh \ 



(81) 



(82) 



(83) 



The fc-dependence of the energy density results of the integration of j^,P{B i ){ r l^iP) with respect to p. It is not 
possible to extend to infinite values of p when we perform the integration, this is because there is a certain scale po 
for which the modes don't have a quantum limit. This can be seen from the mode equation p5|) . where the frequency 
is time dependent 



\+V 2 
if 2 



k 



(84) 



This frequency is always positive, so the solutions oscillate, but depending on the relative magnitudes between p 
and the physical wavenumber k p f lys = ke~ H( ' f , with respect to which of them oscillates. 

In the previous asymptotic limits we have supposed that for modes in the ultraviolet sector, and at the beginning 

of inflation t = U, we have k p h ys {ti) 3> Ho \J \+ p 2 so Ukp{U) — k p h ys (ti)- But, when inflation ends at t = tf, an 

interval of these wavelengths can increase such that k p h y s(tf) <C p and then uJk P (tf) ~ yl /4 + p 2 . This means that 
the actual scales of the observed universe X p hy S will imply a cut po in the integral. The rest of the modes with p > po 
do not have a quantum origin at the start of inflation. 

The connection of inflationary scales to the actual large scales is something uncertain. It is known, approximately, 
the minimum number of e-folds to solve the problems of planarity of the Big-Bang, but this value depends strongly on 
the theoretical model. Furthermore, during reheating there exists an important re-scaling of the cosmological lengths 

M. 



A physical scale today at t = tu, X p h ys corresponds to one at the start of inflation by X phys 



Khys%- In our case 

we are dealing with a inflationary de Sitter expansion, where ^- = e w , such that the scales at the end of inflation 
relate to the actual ones as 



R 



Pinf 



n h radP ; 



h 

cri 



1/4 



(85) 
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N 


H (m p i) 


)& hya (Mpc) 


po «C 


50 


10" 9 


10 4 


— 


60 


10" 9 


10 4 


45 


65 


10~ 9 


10 4 


6.7 x 10 3 


50 


10" 6 


10 4 


— 


60 


10~ 6 


10 4 


1.34 


65 


10" 6 


10 4 


2.1 x 10 2 


50 


1Q -9 


10 2 


— 


60 


10" 9 


10 2 


4.5 x 10 3 


65 


10" 9 


10 2 


6.7 x 10 5 


50 


10" 6 


10 2 




60 


10" 6 


10 2 


1.4 x 10 2 


65 


1Q -8 


10 2 


2.1 x 10 5 



TABLE I: we show some relevant examples on how the parameters N, Hq y Xp hya restrict po- 



the Hubble parameter today is Hh — 2.1h x 10~ 42 GeV, the cosmological radiation density parameter today is 
^rad — h~ 2 4.3 x 10 -5 [32], and the parameter R that describes a phase of reheating [3l| is in units of m p i. It 
doesn't exist in the literature a model in the framework of the Induced Matter theories that studies a realistic 
transition between an inflationary period induced from a 5D vacua and a localized brane dominated by radiation. We 
then assume that this transition is instantaneous. In our model we obtain a vacuum energy density induced in the 
hypersurface 



Pinf 



3 3H 2 



The critical energy density today is 



k 5 % ^5 



8tt 



The parameter R can be expressed as (for instantaneous reheating) 



obtaining the following relation 



lni?= - In pi n f, 

4 



Oh 

a f 



Pinf 



adPc 



1/4 



(86) 



If we use the condition k p hy S {ti) ^ Ho 
with quantum initial conditions 



p 2 we can estimate a cut for the values of p, so we are keeping modes 



Po < 



\ 



„2N 



(87) 



There are three uncertainties in the last relation, the scale of inflation Ho/m p i , the amount of inflation N and the 
scales of coherence of the cosmological magnetic field. Slight variations in this parameters affect significatively the 
cut pq. In tableUwe give some examples on how these parameters N, Hq y \ r p l hys restrict pq. 

Using some numerical calculations [the reader can see the appendix fD])]. we can compute the magnetic en- 
ergy density at a certain scale k, for the modes that propagate on the extra coordinate ■^pm^iVi^iPo) = 

So° dPTgLPmfabp) with 



—P( Bi )(ri,k,po) 



1 



k 4 



2tt 3 H b 5 a ! 



■In 



,1.303 Pott 



8 Hn a 
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where po > A/tt. From the figure (j3|) we see that for reasonable amounts of inflation, N = 60, all of the relevant scales 
are produced. However, for a lower value of e-folds, N — 55, and with a scale of inflation Ho/m p i ~ 10 -6 there are 
excluded the modes corresponding with scales larger than 100 Mpc. 
For the modes Akmiji) we must compute 



1 

1 

2^2 



.4 



km 



The energy density at a certain scale k, is 



dln(k) 



dk k 5 

T¥ 

dk k 5 \A kr . 
k b 5 a 4 



\*A-krn | 



2tt 3 b 5 , 



(89) 



(90) 



Recalling the amplitude (jSTj). we get 



1 



dln(k) 

This also was obtained in [2l|, but in other context. 



2tt 2 6 5 a' 



-H»)(kri) 



in 



(91) 



B. Energy density due to the magnetic tensor Bij 

The energy density that come from magnetic tensor is given by pm,) = — (^ BiJ ^o) 

P(B is) = \ ih.ir- = ^ / **^2 (| +^ 2 ) ^o 2 IA,| 2 , 

from JkP(B^)(v, k,p Q ) = fo° dp-£^p (B%]) {ri,k,p) we obtain 

d , , , 1 A; 2 iJ( 



— p (By) (??,fc,Po) 



27T 3 6 5 a 3 





^ 1.21P07T k \ 


, Po 8" 




\ 8 ffoaj 


2 7T 2 _ 



(92) 



(93) 



where po > 4/V. 

From the modes *4fc m it is easy to see that 



dln{k) 



'b 5 a 2 



(94) 
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If m — — 1 we also have a contribution from the scalar field <f>k in (l52j) and after considering that F(— 1) = 1/2, we 
obtain 



d , , , s 13F,? fc 2 



P(B w )(»7,fc,m =-!) = — 575-2 . (95) 



dln(*)^ w;w ' ' 8Tr 2 b 5 a 

C. Energy density due to the Electric vector Ei 

The energy density related to the electric field 22, it is 

Calculating -^P( Ei )(v, k ,Po) = Jq° d PdM^P(E t )(v, k ,p) w e get 



d . . . 1 fc^iJo 

with po > 4/7r. 



2.77+^ 



where we identify the energy density stored at a certain scale k as 



dln{k) Air 2 b 5 a 4 
this is the same result derived in I2ll. 



D. Energy density due to the Electric scalar E4 



The energy density corresponding to the scalar electric field is P(e 4 ) — ~ (T^^q), so that 



1 r dk ,0 ' 2 
k 



26 3 7 2tt 2 
/" rifc rfp(g 4 ) 

the energy density at a given scale fc, is 

dp(E 4 ) 1 fc 6 



dZn(fe) (9ir 2 H ) 2 b 5 a 6 ' 



(97) 



The energy density associated to the modes Akm is 



^ } - 1 ^ <7(fl)to) M . (99) 



(100) 



(101) 



this component decays strongly with the expansion. It is important to notice that this contribution only appear 
when we are considering a field that doesn't propagate in the extra dimension and has a exponential decay in w with 
rn = —1. 



E. Spectrum of the fields 



In table [XT] we summarize the spectrums of the 5D electromagnetic fields for the two types of fields studied 
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Power spectrum 


5D 


n 


4D 


m< -i 


m>-| 


p 


n 


27r2f/ 0~ 4 I^feP(B) 




~ 5 






J"(m + l)a; 6+2m 




^( M )a; 4+2M 


4 + 2^ 






~ 3 




(m 


+ l) 2 J"(m + l)x 4+2m 


m 2 .F(-m):E 2 - 2m 


/ x 2 J-(M)a; 2+2 ' 4 


2 + 2^ 




■> / „2 \ 
a; 3 f 2.77 + 


~ 3 






g(m)a; 4 + 2m 


g(-m - l)x 2 ~ 2m 


g(6»)x 4+2e 


4 + 26» 








m = —1 










6 



TABLE II: The spectral indexes of the fields that can propagate in the extra dimension are approximate, since there is a 
logarithmic factor and so it doesn't strictly obey a power law. This expressions for the energy density are very blue tilted, so 
they decay to residual values at the end of inflation. The only solutions that can give significant seeds for galactic magnetic 
fields are does arising from real values of m. In particular, when m = —3 we obtain an invariant spectrum for the vector 
magnetic fields. For this value the vector electric and the tensor magnetic fields are red tilted n = —2. For the fields with 
m — —2, the electric and tensor magnetic fields are scale invariant, but the vector magnetic fields is blue tilted with n_g = 2. 



VI. BACK-REACTION EFFECTS 



In order to make an analysis of the back-reaction effects we shall estimate the energy density due to Bij and the 
electric (vector) field, which could be more problematic because, in both cases their energy densities evolve as (kr))~ 2 . 
In order to back-reaction effects to be negligible at the end of inflation, we shall require that these energy densities 
to be smaller than the critical background energy density: p crit — 3HQm 2 [ /(8n). 



From the condition dl £/ k \ P{B i} ) < Pcrit in the eq. (|94|) . we obtain that 



2(1+M) 



3(g /m pi ) 2 ^ 2(1+Al) 



k 

a J 4J 7 (p)p 2 

If we consider p < — 1, modes for which back-reaction effects are negligible has physical wavelengths of the order 

~^ 2 Hp) 



(102) 



^phys\ B < 



V rripi ) 



(103) 



In order to make an estimation we can consider a scale invariant magnetic spectrum with p = — 2 and a Hubble 
parameter Hq = 1.0 x 10 -5 m p i. In this case we obtain that the 



Winy < 0.204 x 10 5 A H , 



(104) 



where A# = 1/Hq is the wavelength related to the Hubble radius during inflation (we are considering c = 1. 

On the other hand, from the eq. (|99|) . we obtain that the condition to back reaction effects do not be relevant in 
the background evolution of the universe, is (for cases with 8 < —2) 



^phys I Ei 



< 



~2g(ey 


2(2 + 8) 


(-) 


3tt 




\m p ij 



A 



H- 



(105) 



that for a scale invariant magnetic spectrum with /i 

A. 



\phys I e 



= m = —2, for which 9 
< 0.512 x 10 5 A H , 



-3, results to be 



(106) 

which is bigger than the Hubble horizon during inflation. Hence, in both cases, for end Ei the modes with 
wavelengths which are on cosmological scales today (these wavelengths are larger than the Hubble radius during 
inflation) are not problematic in the back-reaction sense. On the other hand, it is easy to see that dpB t /[dln{k)] 
is or the order (Ho/m p i) 2 ~ 10~ 10 times smaller than the background energy density, so that this contribution is 
negligible. Of course, in our estimation we are considering a scale invariant spectrum for the magnetic fields, but 
it can be seen from the figure (JSj) that (for Hq ~ 10~ 5 m p ;), the admissible spectrums for the magnetic fields has a 
spectral index —0.12 < nn < 0.1 very close with an scale invariant one. Our results in this sense are in agreement 
with those obtained in [21(. This topic has been discussed in several articles [25l - [27| . 
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VII. SEED MAGNETIC FIELDS 

We are aimed to compare the present day strength of magnetic fields with the predictions of our theory, we can 
estimate the actual magnetic fields generated by this model, the magnetic density parameter defined at a scale 
L = 2-K/k is 

n B (k) = P(B)it = t h ,k) ^ (iQ7) 

Peri 



The energy density related to the magnetic fields after inflation decay adiabatically with the expansion as a 



-4 



P{B) (t = t h , fc) = ^g(t = t f ,k) (108) 



d In k \a,h 
Using the expression (I91[) for the magnetic density, we obtain 

M / „ \ 2 M / U \ 2 ^+ 4 



^1 = (EsU )(*Z±Y (^] " f-^-V . (109) 

dink 9 \m p ij \p m f J \ a // \a h H h J 



Replacing the values for p cr i, pi n f and (a^/a/) we may write [21| 

~ 2.4 x 10 - 7 (2.28 x 10 - 58 r +2 T(p)R-^h^+ 2 ( -^—] ^ . (110) 
dlnfc \a h H h J 

Also, the magnetic density parameter is 

where we use that 1 Tesla 2 /2 ~ 1.9 x 10 -32 GeV A and 1 Tesla = 10 4 Gauss. We then arrive to an expression for the 
seeds of the cosmological magnetic field 



( 77^— ] =s 3.16 x 10- 7 ^ +2 (2.28 x 10~ 58 ) 1+ ^ F(u)^R- , 
\Gauss J \ah Hh 



k ^ I+2 



(112) 



This is the final expression for magnetic fields generated from the fields Akm- The values of R, without considering 
a particular model of reheating, are constrained with a necessary cut for them to preserve nucleosynthesis: 

Pinf > Pnuc - 10~ 85 mp;. 

Furthermore, if inflation holds at sub-Planckian scales we would need that Hi„f/m p i < 1.3 x 10~ 5 , see [3l|. If 
additionally we assume an instantaneous transition from reheating to radiation [which is a good approximation on 
large (cosmological) scales], and keeping in mind that in this case R A ~ (3/8tt)(Hq m 2 ^, we obtain that 

- 24.7 < In R< -5.71, (113) 
which can be written in terms of inflationary scales 

9 x 10~ 43 < — < 1.3 x 10~ 5 . (114) 
rripi 



The scale invariant spectrum in (|112[) is achieved when 



m 



(■^-) ^8.00xl0- 8 (^). (115) 
V Gauss J m= _ 3 \m pi J 



3 For m > —3 the spectrum of magnetic fields is blue tilted. In particular, for m = —1 is 

( B ) ~ 5.09456 X 1(T 65 ( — - — ) 

\GaussJ m= _j \a h H h J 
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The galactic dynamo scenario impose a lower limit for the seed magnetic fields at the moment of the formation of 
seminal galaxies, B > lO~ 22 Gauss (some models may relax to 10 -30 Gauss). This then implies a constraining for the 
scale of the inflation 

— > 1.25 x 1(T 15 . (116) 
rripi 

These values are reasonable for sub-Planckian models. In particular, for Ho/m p i ~ 10~ 9 , we obtain scale invariant 
cosmological magnetic fields of amplitude 



\Gauss J 



.OOxlO -17 , (117) 



m=— 3 



which is compatible with dynamos mechanisms. 

In the figure ([7]) we plotted scale invariant magnetic fields for different inflation scales Hq / m p i versus wavelengths 
expressed in Mpc. The upper constraints come from nucleosynthesis, Faraday Rotation measures and from the Cosmic 
Microwave Background. Furthermore, we include the lower strengths of magnetic fields needed to ignite the galactic 
dynamos and lower magnetic fields from TeV Blazars observations. Notice that for these spectrums the compatible 
scales are: 4 x 10 -7 < H /m p i < 1.3 x 10~ 5 . 

Very recently it was argued a lower bound B~3x 10~ 16 Gauss on the strength of intergalactic magnetic fields, 
which stems from the un-observation of GeV gamma-ray emission from electromagnetic cascade initiated by tera- 
electron volt gamma-ray in an intergalactic medium from Fermi observations in TeV Blazars [33]. More recently, it 
was argued that these results are in tension with inflationary predictions [34|. For our model this implies that the 
minimal scales to produce fields of coherence bigger than Mpc is Hq ~ 4 x 10 _9 m p ;. If the coherent fields are to be 
generated for all scales, including smaller than Mpc then the constraint improves as L -1 / 2 and the scale increases to 
4 x lCT 7 m p ;. 



VIII. CONCLUSIONS 



Using some ideas of IMT, we have studied the magnetogenesis produced by a 5D photon field in a hypersurface 
which undergoes a period of de Sitter inflation. We deduce very blue tilted magnetic fields for photons that can leak 
outside the 4D hypersurface, with wavenumber pHo- In order to produce significative seeds of the magnetic fields to 
explain actual observations, we need other kind of fields that cannot propagate outside the brane. These fields produce 
a discontinuity in the stress tensor, so in turn their origin should come from localized sources (i.e., static solutions) 
on the extra space-like noncompact dimension. We do not study how this fields came from, but once produced we 
analyze their electromagnetic effects in the 4D effective de Sitter (inflationary) universe. We have found that these 
fields are candidates to produce relevant magnetogenesis. In particular, they include the invariant scale spectrum for 
magnetic fields. 

An important result here obtained is that we identify new physical electromagnetic fields that arise from the 
extension of Inflato-electromagnetic Inflation [28|, [3^, [36|, as an extended Maxwell theory. Such fields are, for 4D 
comoving observes in the hypersurface, a scalar electric field and an antisymmetric magnetic tensor field. This 
formalism was proposed recently with the aim to describe, in an unified manner, electromagnetic, gravitational 
and the inflaton fields in the early inflationary universe, from a 5D vacuum. Other conformal symmetry breaking 
mechanisms have been proposed so far [38]. As in our case, most of these are developed in the Coulomb gauge in 
order to simplify the equations of motion for A u . In our calculations, the scalar electric field remains very constrained 
by the actual model, since it is very blue tilted and exists for one specific value of the parameter m. In contrast, 
the magnetic tensor field, Bij, and the vector electric components are red tilted with respect to the vector magnetic 
field. Therefore, its cosmological values at the end of inflation are very important. The electric field is damped during 
reheating giving its energy to charged particles. The tensor magnetic field should also be extinguished, by transferring 
its energy to charged particles or transforming to the vector magnetic field. The final fate of this unobserved quantity 
depends on a model that should explain the transition from this inflationary scenario to the actual universe. However, 
this issue goes beyond the scope of this paper. 

Finally, our results are in very good agreement with recent observations |33j; if the coherent magnetic fields are 
to be generated for all scales, including smaller than AIpc, magnetic fields with strengths larger than 10~ 15 Gauss 
should be generated during inflation on all scales for values of the Hubble parameter larger than Hq > 4 x 10~ 7 m p i, 
which is very compatible with the accepted values during inflation. In the figure © we plotted different admissible 
spectrums of magnetic fields. Notice that if we take into account the more restrictive constrain that comes from TeV 
Blazars, we obtain spectrums nearly scale invariants. The most blue tilted compatible spectrum is for ns — 0.1 and 
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a scale Hq = 1.3 x 10~ 5 m p i. However, the most red tilted magnetic fields take place for ns — —0.62 and a lower scale 
H ~ 10 — 15 77T, p ; . Therefore, taking into account the observational data the spectral indices for magnetic fields should 
be in the range —0.62 < tib < 0.1 for 10~ 15 m p i < H < 10 -5 m p i. 



Appendix A: Conformal coordinates 



From (j4j with the transformation 



t = i> N, 



R = IpoT, 



4> = 4>, 



(Al) 



we obtain the Ponce de Leon metric [37 [ 



Wo 



dt 2 - e 2t ^°dR 2 



d^ 2 



(A2) 



The foliations ip = V'o yield a 4D de Sitter hypersurface with Hubble parameter H = tp^ 1 . 

The vacuum condition R a t, = is widely satisfied because (|A2I) is Riemann flat. On the other hand, the transport 
of vectors in the space is not trivial due to the existence of non zero connections. In coordinates ((31) the non null 
connections are 

^=^\ T) Q =S), Ti = ^T%=e 2N 8), T% = -^e 2N 8). (A3) 

Instead of working with the previous coordinates we will introduce new conformal like coordinates. First we 
introduce a dimensionless coordinate 



(A4) 



yielding 



dt 2 - e 2Hat / Xo dH 2 



dX 2 



(A5) 



The foliation ip — tpo — H Q 1 corresponds in these coordinates with X = Xq = 1. Now we change to a conformal 
coordinate 



3 = lnX, cQ = — = e :t dX, 

X 



(A6) 



such that the metric now reads 



ds 2 = e 2 * 



3 ° (c 



3o 1 dt 2 _ e 2H e-*H 



dK' 



dtf 
Hi 



(A7) 



we can see that taking a constant foliation Z = Zi = cte, we obtain a 4D de Sitter space with Hubble parameter 
H = H e^ 



ds 2 = dT 2 ~e 2HT dv 2 , 



(A8) 



where we have defined the cosmic time as T — e^ 1 ^°t. Here we can perceive that the metric ((H) and their transfor- 
mations generate dS^-spaces for constant foliations in the extra dimension. 

To continue we define the usual conformal time i] — J a(t)~ 1 dt with a(t) = e H ° e °*. If we define the coordinate 
w = Hq 1 ^, the metric finally reads 



ds 2 = b{wf [a{ri) 2 (drj 2 - dR 2 ) - dw 2 ] , 



(A9) 



where the scale factors b(w) — e H ° w y a(rj) — are dimensionless. The components of the coordinates (i], R, w) 

have spatial dimensions. The foliation ip = ipo corresponds with w — wq — in these coordinates, such that b(wo) = 1- 
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Appendix B: Mode quantization 



In order to obtain the canonical structure of ^'A^fa R, wq) and its canonical momentum, we must calculate the 
effective 4D commutators 



= b 3 a 2 



d 



W^faR^^^faR'.tuo) 



(5 M^(r?, R, w ), 7T4m' fa R, w ) 

= tf( 3 '(R-R')S„^ (1+ra+m,) 
The quantization follows in the usual way, by expanding in plane waves e lk R 

d 3 k 



(27r) 3/2 e 4 s =4 

The creation and annihilation operators yield the relations 



The mode equation is 



^A 4m (r h R,w ) = e mH ° 
utors j 

= <5 (3) (k - \a')8 m m', Km,ak 

d 2 (j) k 



a km e ik - R fem fa + aL e " 4kR( / > fc m ( ? ?) 



a km' a k' 



= 0. 



drj 2 



and the normalization condition over them 



- 2?y — h ft fcm = 0, 

ar] 



such that the only accepted value is to = —1, for this model with the Coulomb Gauge choice. 
For the transversal vector field there is no restriction over the values of m 



(5) A J m (77,R, too), fa)-,, faR',™o) 



= b 3 a 2 



( 5 ^„faR,u, ),|^ 5) A im ,faR>') 



iS mm/ e ( - 1+m+m '^ HoWo S { l ) (R - R' 



The mode expansion is given by 
( 5 )A,- m faR,w,„) = e™ H ° tt ' 



d 3 fc 



where the creation and annihilation operators describe an algebra 

& A (k,TO),&l,(k',TO')l =S^(k-k')5 mm ,S xv , [6 A (k,m),6 A ,(k',m / )] = [& A (k,m),& A ,(k',m') 



(Bl) 
(B2) 

(B3) 

(B4) 
(B5) 



(B6) 



]T £jA (km) [6 A (kTO) e 4k - R A fcm (r,) + bl(km)e^ R Al m (fj)] , (B7) 



0. (B8) 



Appendix C: Some properties of the Bessel functions with imaginary order 

The decomposition of the functions (1401) and (|41[) is complex because of their imaginary order. To apply the 
normalization condition we shall use the asymptotic behavior in the short wavelength limit, which describes the 
modes in the UV sector of the spectrum 

(CI) 



J p {x -> +oo) ~ J — cos (x - |) + 0(x~ 3/2 ), 



TTX 
2 



Y p (x -> +oo) ~ J — sin (x - ^ + C(a;" 3/2 ). 



(C2) 



It is interesting to notice that in this limit the modes do not depend in the order p, so in this limit these functions 
become the usual Bessel functions with zero order 

(C3) 
(C4) 



U\;>{x -+ +oo) ~ \l—e~^ ~ n^(x -> +oo), 
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In the short wavelength limit, deep inside the horizon, the fluctuations travel like plane waves in the Minkowski 
spacetime 

A kp \uv — > -j=er lkr >, ki] — -> -oo. (C5) 

V 2fc 

The Hankel function Hp 1 '' (a;) has the desired asymptotic limit. This means that we can set e?2 = and using the 
normalization relationship (J3SJ) , that in this limiting case is simply W[>lfep(r?)] = i, we arrive to d\ = --j^- The 
solution for all scales finally reads 



(C6) 

When the physical wavelengths are much bigger than the Hubble horizon the Bessel functions of imaginary order has 
the following asymptotic behavior 



J p (x -> 0+) 
where 7 p is a phase [3(| defined through 



/2tanh(^) 

PIT 



COS 



tan7 p 



p7rtanh (^) 

Im[r(l + ip)] 



Pln(|)- 
Pln(|) 



7j> 
"7p 



R C [r(i + »p)] ' 



(C7) 
(C8) 

(C9) 



that complies with 70 = 0. This phase was plotted in 



Appendix D: Numerical calculation of ■ykpi'H, k,po) 



Now we must resolve the integral 
ij-'lAkpiv^dp'- 



Po dp 
V 



tanh cos2 



pln(M) 



'7p 



tanh(f; 



We shall resolve approximately this integral. To do it, we separate it in two parts = + From p = 4/ it 
the modes are very close to their asymptotic form 



(Dl) 



V X \A kp (v)\ 2 dp 



4/77 



dp 



P 



tanh ( — ) cos 2 



i' h > ( y J ' v 



pln(^) 



■7p 



tanh (EE) 



The first integral can be numerically resolved 

10& 




dp 

4/tt P ' 



(D2) 



We see that the function approximates linearly to In (j^j- We get 

r,- l \A kp {rj)\ 2 dp ~ In (^j + 1.303. 



The complete integral is then 



The energy density that come from magnetic tensor is given by P(Bi ) = —{T^ 1 



-'\Ak P (v)\ 2 dp^\n(e^°^k V 



0/ 



p {Bij) = \{B ijB «) =^J dkdp^ Q + /) Hl\A kp \\ 

so the energy at a certain scale k with extra momentum p is 

dp Bti _ 1 fc 2 if 2 Q+p 2 ) |2 
dfcdp 2tt 3 65a 2 1 M ' 

As before we split the integral in two parts and solve numerically until the asymptotic regime is r 

4Gn — 




the linear approximation gives 



jTV 1 (~ +P 2 ) \A kp (v)\ 2 dp^ iln + 1.21. 



The energy density related to the electric field -E 1 , it is 

^) ' 1 



2tru 4 



2| /(/ |2 



This energy stored at a scale A: with extra-momentum p is 

rfp(g.) _ 1 fc 2 2 
dfcdp 47r 3 6 5 a 4 ' fepl ■ 

From the expression (|45l) for the amplitude of the temporal derivatives of the modes we get 



dP(Ei) _ 1 pH k 2 



dkdp 7r 3 b 5 a 3 
The numeric integral to p — A/ it gives us 



tanh 



(f) 



sin 



cos 



7p 



tanh (^) 
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FIG. 4: The phase of the complex function T(l + ip) 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



FIG. 5: The amplitude of the modes vs. p for two different scales during the inflationary period in the hypersurface. We can 
see that the modes collapse to p~ for p >> 2/ir. 
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FIG. 6: Amplitude of the derivatives of the modes for two scales. The asymptotic value is 4p 
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FIG. 7: Scale invariant magnetic fields for different inflation scales Ho/m p i plotted against cosmic scales in Mpc. The upper 
constraints come from nucleosynthesis, Faraday rotation measures and from the Cosmic Microwave Background. We also 
include the lower magnetic fields needed to ignite the galactic dynamos and TeV Blazars. For these spectrums the compatible 
scales are: 4 x 10" r < H /m p i < 1.3 x 10" 5 . 
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FIG. 8: Many admissible spectrums of magnetic fields. If we take into account the more restrictive constrain that comes from 
TeV Blazars we obtain that our spectrums should remain close to scale invariance. The most blue tilted compatible spectrum 
is for riB — 0.1 and a scale Ho = 1.3 x 10 _5 m p ;. Furthermore the most red tilted magnetic fields are for ns — —0.62 and a 



lower scale ~ 10 



m v i 



